Dynamics and drivers of mycorrhizal fungi after glacier retreat.

The development of terrestrial ecosystems depends greatly on plant mutualists such as mycorrhizal fungi. The global retreat of glaciers exposes nutrient-poor substrates in extreme environments and provides a unique opportunity to study early successions of mycorrhizal fungi by assessing their dynamics and drivers. We combined environmental DNA metabarcoding and measurements of local conditions to assess the succession of mycorrhizal communities during soil development in 46 glacier forelands around the globe, testing whether dynamics and drivers differ between mycorrhizal types. Mycorrhizal fungi colonized deglaciated areas very quickly (< 10 yr), with arbuscular mycorrhizal fungi tending to become more diverse through time compared to ectomycorrhizal fungi. Both alpha- and beta-diversity of arbuscular mycorrhizal fungi were significantly related to time since glacier retreat and plant communities, while microclimate and primary productivity were more important for ectomycorrhizal fungi. The richness and composition of mycorrhizal communities were also significantly explained by soil chemistry, highlighting the importance of microhabitat for community dynamics. The acceleration of ice melt and the modifications of microclimate forecasted by climate change scenarios are expected to impact the diversity of mycorrhizal partners. These changes could alter the interactions underlying biotic colonization and belowground-aboveground linkages, with multifaceted impacts on soil development and associated ecological processes.

In silico assessment of 18S rDNA metabarcoding markers for the characterization of nematode communities.

Nematodes are keystone actors of soil, freshwater and marine ecosystems, but the complexity of morphological identification has limited broad-scale monitoring of nematode biodiversity. DNA metabarcoding is increasingly used to assess nematode diversity but requires universal primers with high taxonomic coverage and high taxonomic resolution. Several primers have been proposed for the metabarcoding of nematode diversity, many of which target the 18S rRNA gene. In silico analyses have a great potential to assess key parameters of primers, including taxonomic coverage, resolution and specificity. Based on a recently-available reference database, we tested in silico the performance of fourteen commonly used and one newly optimized primer for nematode metabarcoding. Most primers showed very good coverage, amplifying most of the sequences in the reference database, while four markers showed limited coverage. All primers showed good taxonomic resolution. Resolution was particularly good if the aim was the identification of higher-level taxa, such as genera or families. Overall, species-level resolution was higher for primers amplifying long fragments. None of the primers was highly specific for nematodes as, despite some variation, they all amplified a large number of other eukaryotes. Differences in performance across primers highlight the complexity of the choice of markers appropriate for the metabarcoding of nematodes, which depends on a trade-off between taxonomic resolution and the length of amplified fragments. Our in silico analyses provide new insights for the identification of the most appropriate primers, depending on the study goals and the origin of DNA samples. This represents an essential step to design and optimize metabarcoding studies assessing nematode diversity.

The importance of species addition 'versus' replacement varies over succession in plant communities after glacier retreat.

The mechanisms underlying plant succession remain highly debated. Due to the local scope of most studies, we lack a global quantification of the relative importance of species addition 'versus' replacement. We assessed the role of these processes in the variation (ß-diversity) of plant communities colonizing the forelands of 46 retreating glaciers worldwide, using both environmental DNA and traditional surveys. Our findings indicate that addition and replacement concur in determining community changes in deglaciated sites, but their relative importance varied over time. Taxa addition dominated immediately after glacier retreat, as expected in harsh environments, while replacement became more important for late-successional communities. These changes were aligned with total ß-diversity changes, which were more pronounced between early-successional communities than between late-successional communities (>50 yr since glacier retreat). Despite the complexity of community assembly during plant succession, the observed global pattern suggests a generalized shift from the dominance of facilitation and/or stochastic processes in early-successional communities to a predominance of competition later on.

Local climate modulates the development of soil nematode communities after glacier retreat.

The worldwide retreat of glaciers is causing a faster than ever increase in ice-free areas that are leading to the emergence of new ecosystems. Understanding the dynamics of these environments is critical to predicting the consequences of climate change on mountains and at high latitudes. Climatic differences between regions of the world could modulate the emergence of biodiversity and functionality after glacier retreat, yet global tests of this hypothesis are lacking. Nematodes are the most abundant soil animals, with keystone roles in ecosystem functioning, but the lack of global-scale studies limits our understanding of how the taxonomic and functional diversity of nematodes changes during the colonization of proglacial landscapes. We used environmental DNA metabarcoding to characterize nematode communities of 48 glacier forelands from five continents. We assessed how different facets of biodiversity change with the age of deglaciated terrains and tested the hypothesis that colonization patterns are different across forelands with different climatic conditions. Nematodes colonized ice-free areas almost immediately. Both taxonomic and functional richness quickly increased over time, but the increase in nematode diversity was modulated by climate, so that colonization started earlier in forelands with mild summer temperatures. Colder forelands initially hosted poor communities, but the colonization rate then accelerated, eventually leveling biodiversity differences between climatic regimes in the long term. Immediately after glacier retreat, communities were dominated by colonizer taxa with short generation time and r-ecological strategy but community composition shifted through time, with increased frequency of more persister taxa with K-ecological strategy. These changes mostly occurred through the addition of new traits instead of their replacement during succession. The effects of local climate on nematode colonization led to heterogeneous but predictable patterns around the world that likely affect soil communities and overall ecosystem development.

Differential effects of tree species identity on rhizospheric bacterial and fungal community richness and composition across multiple trace element-contaminated sites.

Soil microbial communities play a key role in plant nutrition and stress tolerance. This is particularly true in sites contaminated by trace metals, which often have low fertility and stressful conditions for woody plants in particular. However, we have limited knowledge of the abiotic and biotic factors affecting the richness and composition of microbial communities inhabiting the rhizosphere of plants in contaminated sites. Using high-throughput amplicon sequencing, we studied the rhizospheric bacterial and fungal community structures of 14 woody plant families planted in three contrasting sites contaminated by metals (Pb, Cd, Zn, Mn, Fe, S). The rhizospheric bacterial communities in the given sites showed no significant difference between the various woody species but did differ significantly between sites. The Proteobacteria phylum was dominant, accounting for over 25 % of the overall relative abundance, followed by Actinobacteria, Bacteroidetes and Gemmatimonadetes. Site was also the main driver of fungal community composition, yet unlike bacteria, tree species identity significantly affected fungal communities. The Betulaceae, Salicaceae and Fagaceae families had a high proportion of Basidiomycota, particularly ectomycorrhizal fungi, and the lowest diversity and richness. The other tree families and the unplanted soil harboured a greater abundance of Ascomycota and Mucoromycota. Consequently, for both bacteria and fungi, the site effect significantly impacted their community richness and composition, while the influence of plants on the richness and composition of rhizospheric microbial communities stayed consistent across sites and was dependent on the microbial kingdom. Finally, we highlighted the importance of considering this contrasting response of plant rhizospheric microbial communities in relation to their host identity, particularly to improve assisted revegetation efforts at contaminated sites.

Metabarcoding data reveal vertical multitaxa variation in topsoil communities during the colonization of deglaciated forelands.

Ice-free areas are expanding worldwide due to dramatic glacier shrinkage and are undergoing rapid colonization by multiple lifeforms, thus representing key environments to study ecosystem development. It has been proposed that the colonization dynamics of deglaciated terrains is different between surface and deep soils but that the heterogeneity between communities inhabiting surface and deep soils decreases through time. Nevertheless, tests of this hypothesis remain scarce, and it is unclear whether patterns are consistent among different taxonomic groups. Here, we used environmental DNA metabarcoding to test whether community diversity and composition of six groups (Eukaryota, Bacteria, Mycota, Collembola, Insecta, and Oligochaeta) differ between the surface (0-5 cm) and deeper (7.5-20 cm) soil at different stages of development and across five Alpine glaciers. Taxonomic diversity increased with time since glacier retreat and with soil evolution. The pattern was consistent across groups and soil depths. For Eukaryota and Mycota, alpha-diversity was highest at the surface. Time since glacier retreat explained more variation of community composition than depth. Beta-diversity between surface and deep layers decreased with time since glacier retreat, supporting the hypothesis that the first 20 cm of soil tends to homogenize through time. Several molecular operational taxonomic units of bacteria and fungi were significant indicators of specific depths and/or soil development stages, confirming the strong functional variation of microbial communities through time and depth. The complexity of community patterns highlights the importance of integrating information from multiple taxonomic groups to unravel community variation in response to ongoing global changes.

Optimal sequence similarity thresholds for clustering of molecular operational taxonomic units in DNA metabarcoding studies.

Clustering approaches are pivotal to handle the many sequence variants obtained in DNA metabarcoding data sets, and therefore they have become a key step of metabarcoding analysis pipelines. Clustering often relies on a sequence similarity threshold to gather sequences into molecular operational taxonomic units (MOTUs), each of which ideally represents a homogeneous taxonomic entity (e.g., a species or a genus). However, the choice of the clustering threshold is rarely justified, and its impact on MOTU over-splitting or over-merging even less tested. Here, we evaluated clustering threshold values for several metabarcoding markers under different criteria: limitation of MOTU over-merging, limitation of MOTU over-splitting, and trade-off between over-merging and over-splitting. We extracted sequences from a public database for nine markers, ranging from generalist markers targeting Bacteria or Eukaryota, to more specific markers targeting a class or a subclass (e.g., Insecta, Oligochaeta). Based on the distributions of pairwise sequence similarities within species and within genera, and on the rates of over-splitting and over-merging across different clustering thresholds, we were able to propose threshold values minimizing the risk of over-splitting, that of over-merging, or offering a trade-off between the two risks. For generalist markers, high similarity thresholds (0.96-0.99) are generally appropriate, while more specific markers require lower values (0.85-0.96). These results do not support the use of a fixed clustering threshold. Instead, we advocate careful examination of the most appropriate threshold based on the research objectives, the potential costs of over-splitting and over-merging, and the features of the studied markers.

The generality of cryptic dietary niche differences in diverse large-herbivore assemblages.

Ecological niche differences are necessary for stable species coexistence but are often difficult to discern. Models of dietary niche differentiation in large mammalian herbivores invoke the quality, quantity, and spatiotemporal distribution of plant tissues and growth forms but are agnostic toward food plant species identity. Empirical support for these models is variable, suggesting that additional mechanisms of resource partitioning may be important in sustaining large-herbivore diversity in African savannas. We used DNA metabarcoding to conduct a taxonomically explicit analysis of large-herbivore diets across southeastern Africa, analyzing ∼4,000 fecal samples of 30 species from 10 sites in seven countries over 6 y. We detected 893 food plant taxa from 124 families, but just two families-grasses and legumes-accounted for the majority of herbivore diets. Nonetheless, herbivore species almost invariably partitioned food plant taxa; diet composition differed significantly in 97% of pairwise comparisons between sympatric species, and dissimilarity was pronounced even between the strictest grazers (grass eaters), strictest browsers (nongrass eaters), and closest relatives at each site. Niche differentiation was weakest in an ecosystem recovering from catastrophic defaunation, indicating that food plant partitioning is driven by species interactions, and was stronger at low rainfall, as expected if interspecific competition is a predominant driver. Diets differed more between browsers than grazers, which predictably shaped community organization: Grazer-dominated trophic networks had higher nestedness and lower modularity. That dietary differentiation is structured along taxonomic lines complements prior work on how herbivores partition plant parts and patches and suggests that common mechanisms govern herbivore coexistence and community assembly in savannas.

Environmental DNA of insects and springtails from caves reveals complex processes of eDNA transfer in soils.

Subterranean environments host a substantial amount of biodiversity, however assessing the distribution of species living underground is still extremely challenging. Environmental DNA (eDNA) metabarcoding is a powerful tool to estimate biodiversity in poorly known environments and has excellent performance for soil organisms. Here, we tested 1) whether eDNA metabarcoding from cave soils/sediments allows to successfully detect springtails (Hexapoda: Collembola) and insects (Hexapoda: Insecta); 2) whether eDNA mostly represents autochthonous (cave-dwelling) organisms or it also incorporates information from species living in surface environments; 3) whether eDNA detection probability changes across taxa with different ecology. Environmental DNA metabarcoding analyses detected a large number of Molecular Operational Taxonomic Units (MOTUs) for both insects and springtails. For springtails, detection probability was high, with a substantial proportion of hypogean species, suggesting that eDNA provides good information on the distribution of these organisms in caves. Conversely, for insects most of MOTUs represented taxa living outside caves, and the majority of them represented taxa/organisms living in freshwater environments (Ephemeroptera, Plecoptera and Trichoptera). The eDNA of freshwater insects was particularly abundant in deep sectors of caves, far from the entrance. Furthermore, average detection probability of insects was significantly lower than the one of springtails. This suggests that cave soils/sediments act as "conveyer belts of biodiversity information", possibly because percolating water lead to the accumulation of eDNA of organisms living in nearby areas. Cave soils hold a complex mix of autochthonous and allochthonous eDNA. eDNA provided unprecedented information on the understudied subterranean cave organisms; analyses of detection probability and occupancy can help teasing apart local eDNA from the eDNA representing spatially-integrated biodiversity for whole landscape.

Assessing changes in stream macroinvertebrate communities across ecological gradients using morphological versus DNA metabarcoding approaches.

Freshwater macroinvertebrates provide valuable indicators for biomonitoring ecosystem change in relation to natural and anthropogenic drivers. DNA metabarcoding is an efficient approach for estimating such indicators, but its results may differ from morphotaxonomic approaches traditionally used in biomonitoring. Here we test the hypothesis that despite differences in the number and identity of taxa recorded, both approaches may retrieve comparable patterns of community change, and detect similar ecological gradients influencing such changes. We compared results obtained with morphological identification at family level of macroinvertebrates collected at 80 streams under a Water Framework Directive biomonitoring program in Portugal, with results obtained with metabarcoding from the ethanol preserving the bulk samples, using either single (COI-M19BR2, 16S-Inse01, 18S-Euka02) or multiple markers. Metabarcoding recorded less families and different communities compared to morphotaxonomy, but community sensitivities to disturbance estimated with the IASPT index were more similar across approaches. Spatial variation in local community metrics and the factors influencing such variation were significantly correlated between morphotaxonomy and metabarcoding. After reducing random noise in the dissimilarity matrices, the spatial variation in community composition was also significantly correlated across methods. A dominant gradient of community change was consistently retrieved, and all methods identified a largely similar set of anthropogenic stressors strongly influencing such gradient. Overall, results confirm our initial hypothesis, suggesting that morphotaxonomy and metabarcoding can estimate consistent spatial patterns of community variation and their main drivers. These results are encouraging for macroinvertebrate biomonitoring using metabarcoding approaches, suggesting that they can be intercalibrated with morphotaxonomic approaches to recover equivalent spatial and temporal gradients of ecological change.

Double-digest RAD-sequencing: do pre- and post-sequencing protocol parameters impact biological results?

Next-generation sequencing technologies have opened a new era of research in population genetics. Following these new sequencing opportunities, the use of restriction enzyme-based genotyping techniques, such as restriction site-associated DNA sequencing (RAD-seq) or double-digest RAD-sequencing (ddRAD-seq), has dramatically increased in the last decade. From DNA sampling to SNP calling, the laboratory and bioinformatic parameters of enzyme-based techniques have been investigated in the literature. However, the impact of those parameters on downstream analyses and biological results remains less documented. In this study, we investigated the effects of sevral pre- and post-sequencing settings on ddRAD-seq results for two biological systems: a complex of butterfly species (Coenonympha sp.) and several populations of common beech (Fagus sylvatica). Our results suggest that pre-sequencing parameters (i.e., DNA quantity, number of PCR cycles during library preparation) have a significant impact on the number of recovered reads and SNPs, on the number of unique alleles and on individual heterozygosity. In the same way, we found that post-sequencing settings (i.e., clustering and minimum coverage thresholds) influenced loci reconstruction (e.g., number of loci, mean coverage) and SNP calling (e.g., number of SNPs; heterozygosity) but had only a marginal impact on downstream analyses (e.g., measure of genetic differentiation, estimation of individual admixture, and demographic inferences). In addition, replication analyses confirmed the reproducibility of the ddRAD-seq procedure. Overall, this study assesses the degree of sensitivity of ddRAD-seq data to pre- and post-sequencing protocols, and illustrates its robustness when studying population genetics.

Effects of soil preservation for biodiversity monitoring using environmental DNA.

Environmental DNA (eDNA) metabarcoding is becoming a key tool for biodiversity monitoring over large geographical or taxonomic scales and for elusive taxa such as soil organisms. Increasing sample sizes and interest in remote or extreme areas often require the preservation of soil samples and thus deviations from optimal standardized protocols. However, we still ignore the impact of different methods of soil sample preservation on the results of metabarcoding studies and there is no guideline for best practices so far. Here, we assessed the impact of four methods of soil sample preservation that can be conveniently used also in metabarcoding studies targeting remote or difficult to access areas. Tested methods include: preservation at room temperature for 6 hr, preservation at 4°C for 3 days, desiccation immediately after sampling and preservation for 21 days, and desiccation after 6 hr at room temperature and preservation for 21 days. For each preservation method, we benchmarked resulting estimates of taxon diversity and community composition of three different taxonomic groups (bacteria, fungi and eukaryotes) in three different habitats (forest, river bank and grassland) against results obtained under ideal conditions (i.e., extraction of eDNA immediately after sampling). Overall, the different preservation methods only marginally impaired results and only under certain conditions. When rare taxa were considered, we detected small but significant changes in molecular operational taxonomic units (MOTU) richness of bacteria, fungi and eukaryotes across treatments, but MOTU richness was similar across preservation methods if rare taxa were not considered. All the approaches were able to identify differences in community structure among habitats, and the communities retrieved using the different preservation conditions were extremely similar. We propose guidelines on the selection of the optimal soil sample preservation conditions for metabarcoding studies, depending on the practical constraints, costs and ultimate research goals.

Comparison of markers for the monitoring of freshwater benthic biodiversity through DNA metabarcoding.

Metabarcoding of bulk or environmental DNA has great potential for biomonitoring of freshwater environments. However, successful application of metabarcoding to biodiversity monitoring requires universal primers with high taxonomic coverage that amplify highly variable, short metabarcodes with high taxonomic resolution. Moreover, reliable and extensive reference databases are essential to match the outcome of metabarcoding analyses with available taxonomy and biomonitoring indices. Benthic invertebrates, particularly insects, are key taxa for freshwater bioassessment. Nevertheless, few studies have so far assessed markers for metabarcoding of freshwater macrobenthos. Here we combined in silico and laboratory analyses to test the performance of different markers amplifying regions in the 18S rDNA (Euka02), 16S rDNA (Inse01) and COI (BF1_BR2-COI) genes, and developed an extensive database of benthic macroinvertebrates of France and Europe, with a particular focus on key insect orders (Ephemeroptera, Plecoptera and Trichoptera). Analyses on 1,514 individuals representing different taxa of benthic macroinvertebrates showed very different amplification success across primer combinations. The Euka02 marker showed the highest universality, while the Inse01 marker showed excellent performance for the amplification of insects. BF1_BR2-COI showed the highest resolution, while the resolution of Euka02 was often limited. By combining our data with GenBank information, we developed a curated database including sequences representing 822 genera. The heterogeneous performance of the different primers highlights the complexity in identifying the best markers, and advocates for the integration of multiple metabarcodes for a more comprehensive and accurate understanding of ecological impacts on freshwater biodiversity.

Modelling technical and biological biases in macroinvertebrate community assessment from bulk preservative using multiple metabarcoding markers.

DNA metabarcoding from the ethanol used to store macroinvertebrate bulk samples is a convenient methodological option in molecular biodiversity assessment and biomonitoring of aquatic ecosystems, as it preserves specimens and reduces problems associated with sample sorting. However, this method may be affected by errors and biases, which need to be thoroughly quantified before it can be mainstreamed into biomonitoring programmes. Here, we used 80 unsorted macroinvertebrate samples collected in Portugal under a Water Framework Directive monitoring programme, to compare community diversity and taxonomic composition metrics estimated through morphotaxonomy versus metabarcoding from storage ethanol using three markers (COI-M19BR2, 16S-Inse01 and 18S-Euka02) and a multimarker approach. A preliminary in silico analysis showed that the three markers were adequate for the target taxa, with detection failures related primarily to the lack of adequate barcodes in public databases. Metabarcoding of ethanol samples retrieved far less taxa per site (alpha diversity) than morphotaxonomy, albeit with smaller differences for COI-M19BR2 and the multimarker approach, while estimates of taxa turnover (beta diversity) among sites were similar across methods. Using generalized linear mixed models, we found that after controlling for differences in read coverage across samples, the probability of detection of a taxon was positively related to its proportional abundance, and negatively so to the presence of heavily sclerotized exoskeleton (e.g., Coleoptera). Overall, using our experimental protocol with different template dilutions, the COI marker showed the best performance, but we recommend the use of a multimarker approach to detect a wider range of taxa in freshwater macroinvertebrate samples. Further methodological development and optimization efforts are needed to reduce biases associated with body armouring and rarity in some macroinvertebrate taxa.

Morphological vs. DNA metabarcoding approaches for the evaluation of stream ecological status with benthic invertebrates: Testing different combinations of markers and strategies of data filtering.

Macroinvertebrate assemblages are the most common bioindicators used for stream biomonitoring, yet the standard approach exhibits several time-consuming steps, including the sorting and identification of organisms based on morphological criteria. In this study, we examined if DNA metabarcoding could be used as an efficient molecular-based alternative to the morphology-based monitoring of streams using macroinvertebrates. We compared results achieved with the standard morphological identification of organisms sampled in 18 sites located on 15 French wadeable streams to results obtained with the DNA metabarcoding identification of sorted bulk material of the same macroinvertebrate samples, using read numbers (expressed as relative frequencies) as a proxy for abundances. In particular, we evaluated how combining and filtering metabarcoding data obtained from three different markers (COI: BF1-BR2, 18S: Euka02 and 16S: Inse01) could improve the efficiency of bioassessment. In total, 140 taxa were identified based on morphological criteria, and 127 were identified based on DNA metabarcoding using the three markers, with an overlap of 99 taxa. The threshold values used for sequence filtering based on the "best identity" criterion and the number of reads had an effect on the assessment efficiency of data obtained with each marker. Compared to single marker results, combining data from different markers allowed us to improve the match between biotic index values obtained with the bulk DNA versus morphology-based approaches. Both approaches assigned the same ecological quality class to a majority (86%) of the site sampling events, highlighting both the efficiency of metabarcoding as a biomonitoring tool but also the need for further research to improve this efficiency.

Variability of the Atmospheric PM10 Microbiome in Three Climatic Regions of France.

Primary Biogenic Organic Aerosols (PBOA) were recently shown to be produced by only a few types of microorganisms, emitted by the surrounding vegetation in the case of a regionally homogeneous field site. This study presents the first comprehensive description of the structure and main sources of airborne microbial communities associated with temporal trends in Sugar Compounds (SC) concentrations of PM10 in 3 sites under a climatic gradient in France. By combining sugar chemistry and DNA Metabarcoding approaches, we intended to identify PM10-associated microbial communities and their main sources at three sampling-sites in France, under different climates, during the summer of 2018. This study accounted also for the interannual variability in summer airborne microbial community structure (bacteria and fungi only) associated with PM10-SC concentrations during a 2 consecutive years' survey at one site. Our results showed that temporal changes in PM10-SC in the three sites are associated with the abundance of only a few specific taxa of airborne fungi and bacterial. These taxa differ significantly between the 3 climatic regions studied. The microbial communities structure associated with SC concentrations of PM10 during a consecutive 2-year study remained stable in the rural area. Atmospheric concentration levels of PM10-SC species varied significantly between the 3 study sites, but with no clear difference according to site typology (rural vs. urban), suggesting that SC emissions are related to regional rather than local climatic characteristics. The overall microbial beta diversity in PM10 samples is significantly different from that of the main vegetation around the urban sites studied. This indicates that the airborne microorganisms at these urban sites are not solely from the immediate surrounding vegetation, which contrasts with observations at the scale of a regionally homogeneous rural site in 2017. These results improve our understanding of the spatial behavior of tracers of PBOA emission sources, which need to be better characterized to further implement this important mass fraction of Organic Matter (OM) in Chemical Transport models (CTM).

A heritable subset of the core rumen microbiome dictates dairy cow productivity and emissions.

A 1000-cow study across four European countries was undertaken to understand to what extent ruminant microbiomes can be controlled by the host animal and to identify characteristics of the host rumen microbiome axis that determine productivity and methane emissions. A core rumen microbiome, phylogenetically linked and with a preserved hierarchical structure, was identified. A 39-member subset of the core formed hubs in co-occurrence networks linking microbiome structure to host genetics and phenotype (methane emissions, rumen and blood metabolites, and milk production efficiency). These phenotypes can be predicted from the core microbiome using machine learning algorithms. The heritable core microbes, therefore, present primary targets for rumen manipulation toward sustainable and environmentally friendly agriculture.